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Abstract 

In this letter we investigate how the complex rotation and quivering 
motion of an elongated polarized dust grain in the presence of a monochro- 
matic electromagnetic wave can originate dipolar emission with two dis- 
tinct spectral components. 

We present a model for the emission of radiation by elongated polar- 
ized dust grains under the influence of both an external electromagnetic 
wave and a constant background magnetic field. The dust, exhibiting ro- 
tational motion at the external electromagnetic field frequency u>o as well 
as quivering motion at frequency Qo, proportional to the em field ampli- 
tude, will radiate with frequencies that will depend on the external field 
wavelength and amplitude. 

The radiated spectra exibits a frequency around ojq, and sidebands at 
uJo ± f2o and wo ± 2Qo- Since the amplitude and the frequency of the 
background electromagnetic field are independent parameters, this model 
establishes a correlation between different spectral components of galactic 
dipolar emission, which may help to explain the correlation between a 
component of the Galactic microwave emission and the lOOfim thermal 
emission from interstellar dust, that has been recently measured. 

1 Introduction 

Observations on cosmic background radiation have demonstrated the existence 
of a correlation between a component of the Galactic microwave emission and 
the 100 fim thermal emission from interstellar dust, and several models have 
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been proposed to justify the results measured. Leitcht et al. P suggested free- 
free emission from heated gas to be the cause for this correlation but the model 
could not account for the expected Ha emission unless the temperature would 
be in excess of lO^K, as occurs in shock-heated gas from a supernova remnant. 
Drain and Lazarain showed that for the free-free emission mechanism to be able 
to explain the observed microwave excess it would require an energy input of 
at least two orders of magnitude larger than that provided by supernovae and 
proposed an alternative mechanism to explain the excess of microwave radiation 
observed, based on electric dipole emission from rotating dust grains |21, which 
has been able to account for many aspects of the experimental measurements. 

More recently, data received from the Wilkinson Microwave Anisotropy Probe 
(WMAP), has lead to the claim that the model proposed by Drain and Lazarain 
could only account for 5% of the Galactic microwave emission, and, once again 
rekindled the debate about the origin of " anomalous" dust correlated microwave 
emission 1^. Also, Lazarain and Prunet 0] have analyzed the importance of the 
thermal emission of magnetized dust, locally aligned by the Galactic magnetic 
field in the context of CMB contamination by Galactic dust emission, while 
Ponthieu and others [Bj have stressed the importance of the polarization 
effects by measuring a 3-5% polarized dust signal on the Galactic plane. 

In parallel, several studies of levitation and dynamics of charged dust grains 
have been carried out in both space and laboratory environments [S] , to inves- 
tigate the motion of charged dust particles in low temperature dusty plasmas 
discharges and it has been observed that in a dust plasma sheath dust grains 
levitate due to a balance between gravitational and electrostatic forces, present- 
ing a bouncing motion between electrodes and a quivering motion across the 
electric field of the sheath. 

Several processes have, since then, been proposed to explain dust dynamics, 
with special emphasis on rotational excitation and damping of the dust grains, 
which, in both astrophysical and laboratory plasmas, are known to possess a 
rather elongated shape CHIE] exhibiting a non-zero dipole moment. Examples 
include recoil from thermal collisions and/or evaporation ^2|> collisions with 
gas atoms or plasma and plasma drag absorption and emission of radiation 
[Tl) ^[ , random H2 formation ^21 1 CHI and systematic torques ^7] ^2 ■ How- 
ever, the influence of interstellar electromagnetic background radiation on dust 
dynamics has not yet, to our knowledge, been fully addressed. 

Tskhaya and Shukla jJHI have proposed a simple analytical model to describe 
the dynamics of elongated dust grains in the presence of circularly polarized elec- 
tromagnetic waves and have concluded that the grains not only rotate - or spin 
- with the frequency of electromagnetic fields but also exhibits a complex quiv- 
ering motion. Radiative processes from such spinning dusts have not, however, 
been included in the model. 

This paper addresses the emission of radiation from a single rotating or 
spinning dust grain under the influence of an external electromagnetic (em) 
wave and a constant background magnetic field, correlating the external em 
field amplitude as well as it wavelength with the dust-emitted frequency spectra 
and predicting the radiated power due to both rotation and quivering motions 
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of the dust grain. The work presented is organized as follows: starting from the 
same grain equations of motion as derived in ^2] in Section 2, we address the 
emission of dipolar radiation from spinning dust in Section 3, show the results 
of numerical simulations of the emission spectrum of dust grains in section 4 
and present the conclusions in Section 5. 



2 Spinning dust dynamics 

We start by considering the propagation along the z direction of a circularly 
polarized electromagnetic (em) beam in a medium composed by neutral elon- 
gated dust grains, such as a stellar dust cloud, in the presence of a constant 
background magnetic field B = (0, 0, i?o)- 
The electric field can be written as: 

, t) = Eq exp(z'fc .T')"e (t), (1) 

where i?o is the electric field amplitude and ~e{t) — (cos(woOi sin((jJot), 0) is the 
polarization unit vector. 

We will assume that the dust magnetic moment of each grain is along the 
direction of the background magnetic field (z) and neglect the precession motion 
around this axis, meaning that the dust will rotate in the xy plane. The dipole 
moment of the grain is expressed as d = d(cos </>, sin </>, 0) with </> the orientation 
of the grain relative to the x direction. 

The equation of motion describing the rotation of the grain will be |18| : 

4^ = -AEq sin((/) - woO. (2) 

which can be written as: 

— = -l]oSin((?!)- woi), (3) 

where JIq = dE^ / Iz and Iz is the z component of the principal moment of inertia 
of the dust grain. 

Under the influence of the electromagnetic wave Eq^ the grain will rotate with 
the same angular frequency, ojq d], JH, HOI- Besides this spinning motion, 
the grain may also exhibit a quivering motion corresponding to fluctuations 
on the mean spinning motion. It is therefore convenient to write the angle of 
orientation of the grain as = loqI + 50, which is equivalent to considering 
the quivering motion of the grain in a reference frame which rotates with the 
polarization of the field £^o- The equation of motion then becomes: 

(50 = -f72sin(50). (4) 

As demonstrated in reference |18| . equation Q can be integrated and, im- 
posing the initial conditions: 

6<i>{t = O) = (50o, (5) 
SHt = O) = (50o, (6) 
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results in the following equation: 

2^ ^(50^ - cos(^<?!>) = ^ (^d(j)Q^ - cos((50o) = e (7) 

The constant of integration e plays the role of an effective energy and for 
values in the range [—1,1] the variation of Scj) is bounded, resulting in a quivering 
motion of the grain with a frequency roughly equal to SIq frRj . 

In this paper we are interested in situations when e —1 and the quivering 
motion of the grain is approximately harmonic: 

5(t>it) ~ (500 cos(17o^) + sin(r2ot) (8) 

"0 

In the original frame of reference, the dust grains rotate with an angular 

velocity given hy uj — luq + S(j), thus exhibiting a rotation motion with the 
frequency of the external em field ujq , and a quivering motion across the direction 
of the electric field with a frequency ^Iq, proportional to the external em field 
amplitude. The quivering motion can be identified as the fluctuations of the 
mean rotation motion. 



3 Radiation from spinning dusts 

As the grain rotates about the z axis, the grain emits electromagnetic radiation 
according to the following equation 21 : 



,3^ P( ) 
c2 R 



[it X [it X li]] , (9) 



where the retardation effects have been neglected, "r" is the position vector of 
each element of volume of the grain relative to the center of mass of the grain, R 
is the position vector of the center of mass of the grain relative to the observer, 
n is the direction of observation, p{r ) is the charge density and a is the 
acceleration of of each element of volume of the grain. 

Assuming that the center of mass of the grain remains still the velocity of 
each element of volume of the grain is: 

V = X T*, (10) 

with "cj = (0,0, cjQ + S(f>) the angular velocity vector of the grain, hence the 
corresponding acceleration is: 

li=Saxl^ + Tox [Tu X 1^] , (11) 

where 6a = {0,0,S(/)). The first term in equation (|11|1 corresponds to a tangen- 
tial acceleration at , whereas the second term describes a radial or centripetal 
acceleration a^, which are mutual orthogonal. 
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Using simple textbook algebra, we can derive the following equality: 

[n X [n X ~a]] — — [a^.lt)li + — {a^ .~n)~ri . 

Integrating equation© over the volume of the grain and using the fact that 
R » r, we obtain: 



Erad(y,t) 



R 



6(f>{t)0jf[ut]+i^^{t)0^[^r] 



where we have used the following definitions: 

ul EE (cos((/)),sin((/)),0), 
TC = (— sin((/)), cos(0), 0), 
Oiff[~x] = ~x — {~x .lt)li , 

and where d — J d'^rp{^)~r' is the dipole moment of the grain. 

. 2 

Using the fact that u"^ — uoq + 2ujQ6(j) + 5(j) , yields: 
d 



(12) 



(13) 
(14) 
(15) 



Erad{ r ,t) 



R 



S<j,it)0^[ut] 



LU^ + 2LUoS(j)(t) + d(j){t) 



(16) 

Since the CMB has a low intensity we will investigate the case of weak 
field amplitude, corresponding to a small quivering frequencies and flo << ^^o, 
then the dust-radiated electric field can be separated into three sets of spectral 
lines around luq. The first corresponds to the emission of the standard dipolar 
radiation with frequency loq (remember that ut and Ur rotate approximately 
with frequency uq) associated with the rotational motion of the grain: 



-Brad,l(^,i) 



R 



2 AuJnt 



R 



(17) 



The second set of spectral lines is associated with the quivering motion, 
resulting in a spectral broadening of the standard dipole radiation with the 
generation of two sidebands with frequencies loq + f^o (anti-Stokes) and loq — fig 
(Stokes): 



E rad,2{~r' ,t) 



OC 



R 



R 



[rig + 2a;of^o] argcos(e/2)e' 



(18) 



(19) 
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Finally, the third set of spectral lines is also associated with the quivering 
motion, but the spectral broadening produces essentially two sidebands with 
frequencies coo + 2fio (double anti-Stokes) and uq — 2f2o (double Stokes): 



Erad,z{'r,t) ~ 



R 



oc 



R 



[f72]argcos2(£/2)e*(-°±2no)t. 



The instantaneous energy flux is given by the Poyting vector: 

S = —Erad X Brad = -r- \Eradf ^ , 

and the power radiated per unit solid angle can be written as: 



(20) 
(21) 

(22) 



dP 



An 

dp 



R 



- 25(j) 



\Er 



Hit) 



sin2(et) + 



ujI + 2ujo5(l){t)+6(j){t) 



ujI + 2(JoS(p{t) + 5(t>{t) 



COS(0() C0s(6r) \ , 



(23) 



where &t and 6r are respectively the angle between Ut and Ur, and the direction 
of observation r?. 

By separating the radiation into the two spectral components we recover the 
standard result for the power radiated per unit solid angle emitted by a dipole 
rotating with constant angular velocity luq'- 



dPi d^ 4 . 
^ = ^^oSin (6. 



(24) 



and we identify the power radiated per unit solid angle due to the quivering: 



dQ 



d^ 



2ujo5(t){t) 



— 4:U)oS(l)S(l){t) cos(6t) cos(0r) 



and 



dPs 
dfl 



sin2(er). 



(25) 



(26) 



Since the quivering motion is much slower than the rotation of the grain, 
the average radiated power per unit solid angle over a period of rotation of the 
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Figure 1: Scheme of the mean angular emission of a rotating dust graing: the 
emission is larger along the direction of the background magnetic field and null 
along the plane of rotation of the grain. 



grain yields: 

d{Pi) 



dil Attc^ 

d{P2) _ ^ 

d{Pz) d^ 



a-^[l-cos2($)], (27) 
[nt + 4u;lnl][l-cos'{^)], (28) 



dn - ,^r^^o[l-oos^m], (29) 

where the theorem of equipartion of energy, {6(j){t)) = argcos(£/2), has been 
used and $ representing the angle between the direction of observation "n and 
the plane of rotation of the grain, or the direction of the background magnetic 
field. 

The angular emission distribution is represented in Figure 1, showing that 

the emission is larger along the direction of the background magnetic field and 
null along the plane of rotation of the grain, which indicates that each dust grain 
behaves as a small probe sensing the local electromagnetic environment - both 
the frequency and intensity of the electromagnetic field as well as the direction 
of the background magnetic field - and imprinting them in the emitted dipolar 
radiation. 
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4 Simulation 



The analysis of the dynamics of a single dust grain presented in the previous 
section has provided much insight about the dipolar emission spectrum in a dust 
cloud, enabling the identification of two regimes (for JIq >> l^q and Qq << luq), 
as well as the dominant frequencies of the emitted spectrum. A more real model 
can devised by averaging the emission spectrum over an ensemble of identical 
grains with the quivering temperature T: 

/H= / denis) (/(L^,e)), (30) 

where (/(w, e)) is the Gibbs average of the emission spectrum over the ensemble, 
n(e) ~ noexp [— (e + 1)I{11/2KbT^ the Boltzmann distribution for the proba- 
bilities of occupancy of the rotational energies in the large angular momentum 
limit |22) (notice that the quivering angular momentum of the grain is slflo) 
and Kb the Boltzmann constant. Since the grains are quasi-periodic systems, 
the Gibbs average {I{uj,e)) in equation (|30|) can be replaced by a time sam- 
pling according to the Ergodic hypothesis, which states that the evolution of a 
complex classical dynamical system takes it, with equal probability, through all 
states which are accessible from the starting point subject to the constrain of 
energy conservation |23| |24|. 

The emission spectrum of each dust grain has been simulated using a simple 
fourth-order Runge-Kutta integrator to numerically solve the following equation: 

^ = -sin(<50), (31) 

with T = flot is the natural time scale of the process. To sample over differ- 
ent energy configurations, a method based on the Metropolis algorithm |24| . 
has been used, in which a random configuration of initial conditions is gener- 
ated and the effective energy e is calculated. Should the energy change relative 
to the previous configuration e' be negative, the new configuration is auto- 
matically accepted, else, the new configuration is accepted with probability 
exp [{e-e')ml/2KBT]. 

Simulation results are shown in Figure 2 for both regimes, in which the 
quivering motion is assumed to be in equilibrium with the background at tem- 
perature T = 3K. Also, a cutoff of the effective energy at e = 0.5 has been 
imposed, in order to maintain the results close to the domain of validity of the 
model. The emission pertaining to each single grain is shown in grey (which can 
be interpreted as the sampling fluctuations) and the mean radiated spectrum 
in black. In both regimes the thermal sampling leads to the expected broaden- 
ing of all the emission spectral picks, though it is still possible to identify the 
predicted features. 

Since all nonlinear terms of equation 10} have been included in the sim- 
ulations, it is possible to detect the signature, though weak, of higher order 
sidebands in the fluctuation spectrum (grey lines) at Sloq, Auq... for the low- 
intensity regime, corresponding to nonlinear wave mixing of the external em 
wave through the coupling with the grain. 
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Figure 2: Emission spectrum of the grain: The quivering frequency ilo is much 
smaller than the rotation frequency uo- Sidebands at uq ± fio and coq ± 2Qq are 
clearly identifiable. 
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5 Conclusions 



In summary, we have presented a model for the emission of radiation by elon- 
gated polarized dust grains under the influence of both an external electro- 
magnetic (em) wave and a constant background magnetic field. The emission 
spectrum depends on the wavelength of the external em wave but - and most 
importantly - also on the amplitude, or intensity, of the external em field, which 
is, in fact, the parameter that determines the frequencies of much of the com- 
ponents of the emitted radiation. 

The emission spectrum of the grain exhibits a spectral line centered at the 
external em field wavelength, or frequency, luq, and Raman-like sidebands at 
loq ± r^o and ± 2ilo. We have also identified that the optimum direction of 
emission to be parallel to the background magnetic field, rendering dust grains 
exceptionally useful for local probing of magnetic fields near far astrophysical 
objects. 

Numerical simulations indicate that, as expected, thermal averaging of the 
spectrum, and the inclusion of the nonlinear dust dynamics can lead to the 
broadening and smoothing of the emitted spectrum. 

Since the amplitude and the frequency of the background em field are inde- 
pendent parameters, this model allows to predict a correlation between different 
spectral components of galactic dipolar emission, which may explain the cor- 
relation between a component of the Galactic microwave emission and the 100 
/zm thermal emission from interstellar dust. 
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